Anomalous integer quantum Hall effect in AA-stacked bilayer graphene 
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Recent experiments indicate that AA-stacked bilayer graphenes (BLG) could exist. Since the 
energy bands of the AA-stacked BLG are different from both the monolayer and AB-stacked bilayer 
graphenes, different integer quantum Hall effect in the AA-stacked graphene is expected. We have 
therefore calculated the quantized Hall conductivity a xy and also longitudinal conductivity a xx of 
the AA-stacked BLG within the linear response Kubo formalism. Interestingly, we find that the AA- 
stacked BLG could exhibit both conventional insulating behavior (the v = plateau) and chirality 
for |/x| < t, where v is the filling factor (p = a xy h/e 2 ), p, is the chemical potential, and t is the 
interlayer hopping energy, in striking contrast to the monlayer graphene (MLG) and AB-stacked 
BLG. We also find that for \p,\ ^ [{y/n2 + y/nT) / (y/n~2 — y/ni)]t, where ni = 1, 2, 3, • • •, «2 = 2, 3, 4, • • • 



and ri2 > ni, the Hall conductivity is quantized as a xy — ±^j—n, n — 0, 1,2, • • •, if \fi\ < t 

and a xy = ±^-n, n = 1,2,3, •• •, if |/Z| > t. However, if \p\ = [{y^ + ■s/n^) / {sfna - y/nT)]t, 
the v = ±4(ni + n-z)n plateaus are absent, where n = 1,2, 3, • • •, in comparison with the AB- 
stacked BLG within the two-band approximation. We show that in the low-disorder and high- 
magnetic-field regime, o xx — > as long as the Fermi level is not close to a Dirac point, where F 
denotes the Landau level broadening induced by disorder. Furthermore, when a xy is plotted as a 
function of p,, a v = plateau appears across p = and it would disappear if the magnetic field 



B = ■Kf/Nehvf, N = 1,2,3, • 



Finally, the disappearance of the zero-Hall conductivity plateau 



is always accompanied by the occurence of a 8e //i-step at fi = t. 
PACS numbers: 73.43.Cd,71.70.Di,72.80.Vp,73.22.Pr 
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I. INTRODUCTION 

Graphene exhibits many peculiar properties [l[ and has 
greatly intrigued physicists in recent years. Charge carri- 
ers in the monolayer graphene (MLG) possess a linear en- 
ergy dispersion (see Fig. la) @ and are of chiral natureQ 
near each Dirac point. The quasiparticles in the AB- 
stacked bilayer graphene (BLG) are also chiral. However, 
unlike MLG, the energy spectra of the AB-stacked BLG 
are parabolic (Fig. lb). One of the interesting properties 
of graphene is quantum Hall effect (QHE). Indeed, both 
theoretical^ and experimental works0-0] show that in- 
teger quantum Hall effect (IQHE) in MLG is unconven- 
tional. The Hall conductivity in MLG is quantized as 
<j xy = ±^|-(n + \), where n = 0, 1,2, • • •. The factor 
4 comes from the fourfold (spin and valley) degeneracy. 
Furthermore, beacuse the states are shared by electron 
and hole at the zeroth Landau level (LL), the shift of 1/2 
occurs. In contrast, in the AB-stacked BLG, within the 
two-band parabolic approximation, the Hall conductivity 
was shown to be a xy = ±^-n, where n = 1,2,3, • ■ •[!]■ 
The zeroth and first Landau levels are degenerate and 
hence the first quantum Hall plateau appears at 4e 2 //i 
instead of 2e 2 /h. This phenomenon has been observed 
experimentally QHE of the AB-stacked BLG was also 
studied based on a four-band Hamiltoian in Ref. 11011 . 



'Electronic address: gyguo@phys.ntu.edu.tw 



Although AB stacking is predicted to be energetically 
favored over AA stacking in ab initio density functional 
theory (DFT) calculations, the energy difference of about 
0.02 eV/cell is small pi El]. Moreover, Lauffer et al. 
found that scanning tunneling microscopy (STM) images 
of BLG resemble that of MLG, and hence they regarded 
it as a consequence of the BLG confiquration being close 
to AA stackingil3|. Moreover, Liu et al. reported that 
in their high-resolution transmission electron microscope 
(HR-TEM) experiments a high proportion of thermally 
treated samples are AA-stacked BLG [l4[. These findings 
indicate the possibilty of fabrication of the AA-stacked 
BLG. Since the energy bands of the AA-stacked BLG (see 
Fig. lc) are different from both the AB-stacked BLG and 
monolayer graphene, the quantum Hall effect in the AA- 
stacked is expected to be quite different from that in the 
latter two systems. 

We have therefore carried out a theoretical study of 
IQHE as well as the longitudinal conductivity a xx in the 
AA-stacked bilayer graphene using the Kubo formalism. 
In this paper, we present a general analytical form of 
the Hall conductivity [a xy {jl, B)] as a function of both 
chemical potential (/2) and magnetic field (B) of the AA- 
stacked BLG. Our presentation will be divided into two 
parts: i) the variation of the a xy vs. 1/B curve with 
some fixed /2's and ii) the effect of the magnetic field on 
the o~ xy vs. [2 curve. Our main findings follows. 
Firstly, for \p,\ ^ [(^/n^ + y/hT) / (y/n^ - y/h~T)]t, where t is 
the interlayer hopping energy, m and n-x are any integers 
larger than 1 and 2, respectively, and n\ < n%, the Hall 
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FIG. 1: (color online) The energy bands of (a) mono- 
layer graphene, (b) AB-stacked bilayer graphene, and (c)AA- 
stacked bilayer graphene. 



conductivity is quantized as 

4e 2 

a xy = ±—n, n = 0, 1, 2, • • •, if < t 
4e 2 

a xy = ±—n, n = 1, 2, 3, • • ■, if \fi\ > t. (1) 



However, if = [(y 7 ^ - 
conductivity is given by 



ni)/(y/ri2 — y/rii)]t, the Hall 



, 4e 2 
±— n, 
h 



4e 2 
~h 



excluding ±——(n 1 +n2)n, n— 1,2,3, 



(2) 

Secondly, in the low-disorder and high-magnetic-field 
regime [r 

2 M)[(/i 2 



and h 4 ujf > 



i 2 )r 2 /^ c 4 )] 



(/i + Vt)% <7 XX RJ 

i 2 )r 2 /(/2 2 - i 2 ) 2 + 2T 2 /h 2 oj 2 + 5(/2 2 + 
► if the Fermi level is not close to a 
Dirac point, where Huj c is the cyclotron energy. That is 
to say, if the magnetic field is high enough, the applied 
electric field cannot drive any current for |/2| < t, while 
the current is perpendicular to external electric field for 
\fi\ > t. Thirdly, we find that the a xy = (the fill- 
ing factor v = 0) plateau across p, = would disap- 
pear when B = irt 2 /Nehvp, N = 1, 2, 3, • • •, and that a 
8e 2 //i-step at p, = t occurs while a zero- Hall conductivity 
plateau disappears. Interestingly, unlike the monolayer 
and AB-stacked bilayer graphenes, the AA-stacked bi- 
layer graphene could display a unusual v = plateau 
even though it contains chiral quasiparticles. We argue 
that the occurrence of the v = plateau is due to the 
shift of level anomalies by the interlayer hopping energy 
t. 



II. THEORETICAL MODEL AND 
ANALYTICAL CALCULATION 

A. Model Hamiltonian and Landau levels 

We first derive an effective four-band Hamiltonian near 
each Dirac point for the AA-stacked bilayer graphene in 
tight-binding aprroximation via k • p expansion We 
then diagonalizc this Hamiltonian to obtain the energy 
bands of the AA-stacked BLG. In Fig. 1, we display the 



energy bands of the MLG, AB-stacked, and AA-stacked 
BLG together. We find that the energy bands of the 
AA-stacked BLG are just two copies of the MLG band 
structure shifted up and down by t, respectively. Hence, 
E = ±t are the Dirac points for the AA-stacked BLG. 
When a magnetic field is applied, we should replace the 
momentum operator p with p + eA/c, where the exter- 
nal magnetic field B = V x A and -c is the charge of an 
electron. The magnetic field is applied along the positive 
z-axis (i.e. out of plane) and hence the vector poten- 
tial can be written as A = (—By ,0,0) in the Landau 
gauge. Therefore, the effective four-band Hamiltonian in 
the presence of the magnetic field is given by 



H+ = 



Vf{<J x TT x ± <JyTTy) 
-tl 



-tl 



ihd T 



eBy/c, 



TTy = —ihdy. 



(3) 



Here ± label the two valleys of the band structure at K 
and K' , respectively, vf denotes the Fermi velocity. In 
the Laudau gauge, we can substitute the eigenfunction 



v> = 



'4>(y) of the Hamiltonian into the schrodingcr 



equation Hip = 

(My),My)) T : 

column vectors. 
a ± = a x ±ia y , £ 



Etp. Here 4>(y) can be written as 
where 4>i and 4>2 are two-component 
Then, we make the transformations: 
= y/e B - e B k and OT = (£± d c )/y/2 , 

Finally, for 



where the magnetic length £b = y/hc/\eB\. 
the K valley, the schrodinger equation reads 



V2£ E 



(0-a+ H 
-tl 



O+a' 



—hyp 



-tl 

(0-(T+ + 0+a-) 



E 



01 

<P2 



(4) 



Since T satisfy the commutation relation: [0~ , + } = 
1, are the annihilation and creation operators of one- 
dimensional (1-D) simple harmonic oscillator (SHO), re- 
spectively. Similarly, er^ are the raising and lowering 
operators of pesudospin angular momentum. Obviously, 
the eigenstates of 0~a + + + a^ are (\N - 1), ±\N)) T if 
N > 1 and (0, |0)) T if N = 0, where \N) arc the eigen- 
states of the 1-D SHO. All non-zero vectors are eigenvec- 
tors of —tl. Therefore, we can infer that the eigenvalues 
of Eq. (4) are 

Eft = -fiVNhuj c - vt (5) 
with oj c = V%vf/£b and the eigenstates of Eq. (4) are 



Hjv) 




if N = 0, 



if N > 1, 



(G) 
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Here the indicc p = ± and v = ±. Clearly, the LLs of 
the AA-stacked BLG are just two copies of the LLs of 
the MLG shifted up and down by t, respectively. 



B. Linear response calculation 

The conductivity can be calculated using the Kubo 
formula within the linear response theory [7j|. The Kubo 
formula for the DC-conductivity is given by 



lim 

n->o 



Imng(ft + iO) 



(7) 



Here the retarded current-current correlation in the 



Matsubara form reads (To| 
4e 2 



(8) 



E EEE 



(k, p, v\vi\e, p, a) (£, p, a\vj\k, p, v) 



(iU! n 



ED 



where the factor 4 is due to the fourfold (spin and val- 
ley) degeneracy, the velocity operator Vi = [xi H+]/ih, 
and E£ v — El" /K. u> n and v m arc Matsubara fre- 
quencies of fermion and boson, respectively. When the 
chemical potential and disorder scattering are considered, 
Matsubara frequency of fermion has to be corrected as 
iuj n = iu n + p/h + isgn(w„)r / H. T is the Landau level 
broadening due to the presence of disorder. Substituting 
the eigenstates of H + into Eq. (8), we obtain 



where \k t 



2„,2 



le v 



(2E^- 



is defined as 



(9) 



£>1 v=± 



n—oo 

Ei 

n— — oo 
{iVm ~ 



(iUJ n - ^ k 

► —iv m )]. 



Ei v ){iw 



E 



i>"\ 



(10) 



In the DC and clean limit, we let ft + iO — > and set 
L = 0. Then, after evaluating the Matsubara sumsflil]. 
we find that 

„ -2(fi + iO)[/(£T)-/^r)] 

"Xfc l |w m =n+iO->0 ~ 



1 

pk' 



(11) 

where f(E) is the Fermi-Dirac distribution. Further- 
more, we define f(E) = f(E) — 1/2. Then, using Eqs. 
(7), (9) and (11) and considering the variation of direc- 
tion of conductivity with the signs of magnet field and 
carrier charge, we can derive that 



>xy 



4e 2 

-— sgn(e£) 



(12) 



At zero temperature, f(E) = 1 and f(E) = for the 
occupied and unoccupied LLs, respectively. The LLs lo- 
cated at E < p are ocuppied and the others are empty. 
That is to say, f(E) = 1/2 for E < p, while f(E) = -1/2 
for E > p. Therefore, we only need to calculate the num- 
ber of LLs between —p and p to determine the magnitude 
of the Hall conductivity. Moreover, for —\p\ < E < \p\, 
the number of the up-shifted LLs (y = — ) is equal to 
that of the down-shifted LLs (y = +). Hence, we find 
that the zero-temperature Hall conductivity is given by 



4e" 



• xy 



sgn(p)sgn(eB) [9(\p\ + t)d(\p\ - t) 



e>i n=± 



(13) 



Eq. (13) indicates that the Hall conductivity would sim- 
ply change its sign as p — > — p and that the Hall conduc- 
tivity is equal to Ae 2 /h times the number of up-shifted 
LLs between —\p\ and \p\. Eq. (13) can also be written 
in the form 



4e 2 

Txy = — — sgn(//)sgn(e5)x 



p\-t-j2hv 2 F \eB\/c) 



C (\p\-ty 



9{j2hv 2 F \eB\/c-t+\p\ 



2hv 2 F \eB\ 
'c{\p\-tf 



2hv 2 F \eB\ 



(14) 



\c(\p\- 


f t) 2 ' 




\c(\p\- 


-t) 2 ] 


2hv 2 F 


eB\_ 




2hv 2 F 


eB\\ 


K\P\ + 


t)9{\p\ - 


t)}. 





Here [x] means the integer part of x. The last term is 
the contribution of level anomalies (N = 0)Q. 

We also calculate the longitudinal conductivity (cr X x) 
via the Kubo formula. The longitudinal conductivity 
have to be evaluated under the effect of disorder scat- 
tering (T #0). Using Cauchy's integral theorem as in 
Refs. [mES, we can obtain 



nh 



E E E Im Ci(£) Im 5f (E) + 2j2^gr(E)lm~g^(E) 

(15) 



1 1 ,u 

=± 



g' e = l/(E/h- J5f + iT/h). Applying the tech- 
niques of partial-fraction decomposition similar to that 
used in Ref. [l8j and after some cumbersome algebra, we 
finally find that the zero-temperature longitudinal con- 
ductivity can be written in terms of the digamma func- 



tion (^0)0, 

E 



4e^ 
irh 



■-PQ 



± 

(ft + Vtf 



(p + vtf 



r 2 2Y(n 



[2(/2 + i4) 3 r + 2(/2 + j/<)r 3 ] 



+ 1 



vt) 



2T(p + id) 

1 Wj 2 



jl = 0.10eV 
(i = 0.15 eV 
jx = 0.25 eV 
\l = 0.30 eV 



p 



[h 6 uj 6 c + rf^ip + vtf - nh 2 Lu 2 c (p + vtf + 8(/2 + vtf] r 2 

q 

[SV + 4ft 2 w 2 (/2 + vtf + 16(m + vtf] r 4 



8{p + vtfr 6 
q 



vtf + h 4 uii] r 2 + h 2 uj 2 T 4 



[{p. + hw c + vtf + r 2 ] [(p - hu) c + vtf + r 2 ] 
1 



(P + vtf + r 2 

Here, 



(16) 



t4 4 



P 

q = {[(fi 



- i6r 2 (/i 

vt) 



vtf 



h 2 uj 2 c - r 2 ] 2 + 4T 2 {p + vtf} x p. 



(17) 



Unlike Ref. 12CJ, we do not adopt low-magnetic field 
approximation here and therefore Eq. (16) is general 
and suitable for all values of the applied magnetic field. 
However, when we adopt the low-disorder and high- 
magnetic-field limit, we first consider T as being small 
and keep the terms in the order of T 2 . Then, we let 
(fi + vt)/(h 4 uj 4 ) —> 0, Eq. (16) could be simplified as 



4e^ 
■nh 



(P 2 + t 2 )r 2 
(p 2 - 1 2 ) 2 



r 2 



h 2 L 



. p 2 + t 2 T 2 
h 2 tu 2 h 2 uj 2 



(18) 

The first term is independent of B. The second and third 
terms are proportional to 1/B and 1/B 2 , respectively. 



III. DEPENDENCE OF CONDUCTIVITY ON 
MAGNETIC FIELD AND CHEMICAL 
POTENTIAL 

Since both the magnetic field and chemical potential 
could be tuned experimentally, we display the calculated 
conductivity as a function of 1/B and p in this Sec. Here 
we use the interlayer hopping energy t = 0.2 eV for the 
AA-stacking and t = 0.4 cV for the AB-stacking as deter- 
mined by the ab initio DFT calculations within the local 
density approximation (LDA) [2l| . Moreover, because the 
quantized values of the Hall conductivity is independent 
of the presence of disorder scattering [22j, we show only 
the Hall conductivity in the clean limit (i.e., T = 0) and 
analyze the Hall plateaus qualitatively. 
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FIG. 2: (color online) (a) The quantized Hall conductivity 
a xy of the AA-stacked bilayer graphene as a function of 1/B 
for several values of chemical potential p. (b) The quantized 
Hall conductivity a xy of both AA-stacked and AB-stacked bi- 
layer graphenes as a function of 1/B. The interlayer hopping 
energy t used is 0.2 eV for the AA-stacking and 0.4 eV for the 
AB-stacking, and vf = 1.0 x 10 6 m/s. The Hall conductivity 
<j xy for the AB-stacking was obtained by using Eq. 15 from 
Ref. 113. 



A. Conductivities vs inverse of magnetic field 

In Fig. 2, the Hall conductivity is plotted as a func- 
tion of 1/B and we should discuss the effect of chem- 
ical potential on the a xy vs. 1/B curve. From Fig. 

2(a), we see that a xy — ±^n, with n = 0, 1,2,- • •, 
for \p\ < t, and n = 1,2,3, • • •, for \p\ > t, excluding 
|A*| = t(V^2 + \/n{)l (V"2 - ^/ni)]t, where m = 1, 2, 3, ■ ■ ■ 
and ri2 = 2, 3, 4, It is clear from either Eq. (13) or Eq. 
(14) that <J X y{—p) = —&xy(P)i an d hence we did not show 
any curves for p < in Fig. 2(a). Interestingly, in con- 
strast to the MLG and AB-stacked BLG, the AA-stacked 
BLG displays the pronounced v = plateau for \p\ < t, 
where the filling factor v = a xy h/e 2 . The MLG and AB- 
stacked BLG lack the v = plateau because their level 
anomalies are located at E — 0. The level anomaly of the 
MLG is the zeroth Landau level while those of the AB- 
stacked BLG are the zeroth and first Landau levels. The 
occurrence of level anomalies is a remarkable manifesta- 
tion of the unique property of chiral quasiparticles [8l [23| . 
In the AA-stacked BLG, similarly, level anomalies also 
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two 8e 2 //i-steps, the curve passes through three Ae 2 /h- 
steps and four plateaus. Only the Hall plateaus a xy = 



0.08 



0.04 
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FIG. 3: (color online) The longitudinal (a xx ) and transverse 
(cr X y) conductivities of the AA-stacked bilayer graphene as a 
fuction of 1/B, (a) for \p,\ < t, and (b) for \p,\ > t. The rest 
parameters are the same as in Fig. 2. 



exist in the LL spectrum. However, for the AA-stacked 
BLG, the level anomalies are shifted up and down by 
t, respectively. Thus, these level anomalies are unique in 
the sense that they are always located at the Dirac points 
regardless of the magnitude of B, in contrast with the 
other Landau levels jit l8llloj|. It is worth mentioning that 
the Dirac points often exhibit interesting electronic prop- 
erties, such as electron- hole puddle formation [U [H|, and 
Andreev reflection type transitions. [26| Untill now, some 
transport properties at these Dirac points remain to be 
understood [27[. The level anomaly is one of the interest- 
ing properties of the Dirac points and recently the nature 
of its electronic states (being metallic or insulating) in 
the high magnetic field-low temperature regime is hotly 
debated [28l|. 

Displayed in Fig. 2(b) are the Hall plateaus of the 
AA-stacked BLG for \p\ = [(y/n^+ y/m)/(\/n2- \/ni)}t. 
We find that in addition to the v = plateau, other dif- 
ferences between the AA-stacked and AB-stackcd BLGs 
exist. In particular, when \p\ = [(v/nij + -^ni) / {y/ni — 
y/nT)]t, some 8e 2 //i-steps appear at 1/B ^ 0. In other 
words, in comparison with the other cases of \p\ > t 
[(v / "i + v / "2)/(v / "2 — some plateaus would 

be missing for \p\ = [(^/nj + y/fii) / (y/fi® - y/ni)]t. Fur- 
thermore, these 8e 2 //i-steps appear periodically. Taking 
p = 3t (i.e. ri2 = 4, m = 1), for example, between any 



4c- 



±1, ±2, ±3, ±4, ±6, ±7, ±8, ±9, ±11, ±12, ■ 



•, appear. In other words, the Hall plateaus n = 
±5, ±10, ±15, • • • are absent here. When \p\ > t, the 
quantum Hall effect of the AB-stacked BLG must be 
studied based on a four-band Hamiltonian[lfl| . Based 
on the four-band model[Ioj], for \p\ > t, the AB-stacked 
BLG can also exhibit a 8e 2 //i-step, as shown in Fig. 
2(b). Although the 8e 2 //i-step is not specific to the AA- 
stacked BLG, the periodic appearance of the 8e 2 //i-steps 
has never been seen in the AB-stacked BLG and hence 
is a unique characteristic of the AA-stacked BLG. 

The longitudinal and transverse conductivities as a fuc- 
tion of 1/B are plotted together in Fig. 3. It is seen 
from Fig. 3 that the longitudinal conductivity goes to 
a local minima at the position of the Hall plateaus and 
reaches a local maxima as the steps appear except the 
step at 1/B = 0. The unique v = Hall plateau of 
the AA-stacked bilayer graphene for \p\ < t is especially 
interesting. From Fig. 3(a), we find that a xx falls to 
zero as the v = Hall plateau emerges. As stated previ- 
ously, for \p\ > t, the AA-stacked BLG lacks the v = 
Hall plateau. The first Hall plateau occurs at v = 4 or 
v = —4. This encourages us to investigate the difference 
between the longitudinal conductivities at the v = Hall 
plateau for |/2| < t and at the v = ±4 Hall plateaus 
for \p\ > t. Fig. 3(b) shows that for |/2| > t, the <r xx 
goes to zero as the first Hall plateau occurs. This implies 
that at the high magnetic field, the external electric field 
cannot drive any current for \p,\ < t while the current 
is perpendicular to the external electric field for \p,\ > t. 
Here Jl and t are in the order of 0.1 eV while the order 
of magnitude of T [0(T)] is 0.01 eV. p, + vt are about 
one order of magnitude higher than T (fl + vt ~ 10T). 
Therefore, the condition of low disorder is satisfied. If 
h 4 ^ ^> (p, + vt), Eq. (18) can be applied here. This 
needs Huj c to be larger than 1.78 (p, + vt) and B is esti- 
mated to be at least a few times larger than 10 T. In the 
low-disorder and high-magnctic-ficld regime, we roughly 
estimate from Eq. (18) 



4^ 



(19) 



When 0(\p\ - 1) > IT, a xx - 0.1(4e 2 //i) 0. 

Interestingly, the results shown in Fig. 2 could be 
explained in terms of Fig. 4. Let us account for Fig. 
2(a) first. It is clear from Fig. 4 that for \p\ < t, level 
anomalies are outside the range of — \p\ ~ \p\. Con- 
versely, for \p\ > t, level anomalies arc inside the range of 
— \p\ ~ \p\- In the high magnetic field regime, for \p\ < t, 
no Landau level exists between —\p\ and \p\ and hence 
the AA-stacked BLG displays the conventional insulating 
behaviour (a v = Hall plateau) even though it possesses 
chirality. Such behaviour of the AA-stacked BLG is in 
stark contrast to the MLG and AB-stacked BLG. How- 
ever, for \p\ > t, level anomalies are located between — \p\ 
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FIG. 4: (color) The Landau level spectrum of the AA-stacked 
bilayer graphene (a) with \fi\ < t and (b) with \p,\ > t, re- 
spectively. Here circles mark the positions of Landau lev- 
els. Red, green, and purple circles represent the locations of 
level anomalies as well as other up-shifted and down-shifted 
Landau levels, respectively. Mazarine and blue lines denote 
the up-shifted and down-shifted energy bands, respectively. 
Shaded is the region between —p, and \p,\. 

and 1/2 1 even in the high magnetic field regime and hence 
a v = plateau cannot emerge. 

The Hall plateaus displayed in Fig. 2(b) can be ex- 
plained as follows. A up-shifted LL of -E/T M ' _ and a down- 
shifted LL of E£ ,+ are partners because Et = —E7 tl ~. 
As the up-shifted LL of E^* 1 ' - goes through the |/2|-level, 
the down-shifted LL of E%' + passes through the — \p\- 
level. They always enter the region between |/2| and 
(i.e. the shaded region) together and hence each con- 
tribute 4e 2 //i to the Hall conductivity. Therefore, the 
Hall conductivity is equal to 4e 2 /h times the number of 
cither up-shifted or down-shifted LLs between — 1/2| and 
\p,\. Therefore, in order to form a 8e 2 //i-step, either two 
up-shifted or down-shifted LLs must enter or leave the 
shaded region together. Hence we only need to focus on 
either up-shifted or down-shifted Dirac cones and discuss 
the movement of the LLs located in this cone to explain 
the origin of 8e 2 //i-stcps. Let us consider the up-shifted 
Dirac cone. As the magnetic field decreases gradually, the 
up-shifted LLs would go close to its level anomaly, E = t. 
Therefore, we can infer that for |/2| < t, the LLs above 
the |/2|-level (called the upper LLs) go far away from the 
|/2|-level, while the LLs below the — |/2|-level (called the 
lower LLs) move toward the — 1/2| -level, as shown in Fig. 
4(a). The lower LLs would enter the shaded region one 
by one but the upper LLs can never get into the shaded 
region. However, for |/2| > t, both the upper and lower 
LLs can go close to the shaded region [see Fig. 4(b)], i.e., 
two up-shifted LLs may enter the shaded region together. 
Therefore, the 8e 2 //i-steps can only appear when \p\ > t. 

For brevity, we use indices (k, /i, v) to denote the Lan- 
dau level of Efc v below. Let us label the two LLs which 
enter the shaded region together as the (fci,— ,— ) and 
(*&,+,-) LLs. Then, k Y jk 2 = (|/2| - t) 2 /(|/2| + t) 2 . 



TABLE I: The characteristics of integer quantum Hall effect in 
monolayer (ML) , AB-stacked bilayer ( AB) and AA-stacked bi- 
layer (AA) graphenes as well as conventional two-dimensional 
semiconductor structures (2D). 





ML 


AB 


AA 


2D 


plateau steps (5-) 


4 


4, 8 a 


4, 8 6 


2 


v = plateau 


no 


no 


yes c 


yes 



"The 8^j- plateau step only occurs (aperiodically) in the four band 
model (Ref. [Toj). 

h 2 

The 8^- plateau step appears periodically for [ju| > t only. 
The v = plateau occurs only if \p\ < t. 
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FIG. 5: (color online) The quantized Hall conductivity a xy 
of the AA-stacked bilayer graphene as a function of chemical 
potential ft for several values of magnetic field B. The rest 
parameters are the same as in Fig. 2. 



(711,712) satisfies this condition and ni/712 is an irre- 
ducible fraction. Then, (fci,fc 2 ) = (pni,pn 2 ), where p = 
1,2, 3--, is a set of solutions of fci/fc 2 = t) 2 / (\fi\+t) 2 . 
Between the entries of [{p— l)ni, (p— 1)712] and [pni,pri 2 ] 
LLs, (ni +712 — 2) LLs get into the shaded region sequen- 
tially as 1/B decreases. Hence, (ni +n 2 — 2) 4e 2 //i-steps 
occur between any two 8e 2 //i-steps. The Hall conduc- 
tivity is quantized as a xy = ±H-n with the exception 

of ±^|-(ni + 712)71, where n = 1, 2, 3, • ■ ■, as shown in 
Fig. 2(b). Unlike the AB-stacked BLG and the other 
cases of 1/2 1 > t of the AA-stacked BLG, the Hall con- 
ductivity for |/2| = [(^/na + \f^) / (\fn2 — y/ni)]t lacks 
the v = ±4(rii + n 2 )n plateaus. Furthermore, it is clear 
from Fig. 2(b) that when \p\ > t, the AB-stacked BLG 
can also exhibit a 8e 2 //i-step but the appearance of the 
8e 2 //i-steps is not periodical. The main findings here are 
summarized in Table I. 



B. Hall conductivity vs chemical potential 

Fig. 5 is a plot of a xy versus /2, showing how the Hall 
plateaus are influenced by the magnetic field. It is clear 
from Fig. 4 that unlike the MLG and AB-stacked BLG, 
a v — plateau centered at p, = appears for B = 12 
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T in the AA-stacked BLG. However, when the condition 
of \/2NhvF / (-B — t is reached by tuning the magnetic 
field, the (N,+,—) and (N, — ,+) LLs would be located 
at E = 0. As a result, the v = plateau disappears and a 
8e 2 //i-step at p = forms, like the AB-stacked BLG. In 
other words, the absence of the v = plateau needs the 
magnetic field B = nt 2 /Nehvp, where N = 1,2 • ••. In 
Fig. 5, B = 10.1 T and S = 15.2 T satisfy this condition 
with N = 3 and N = 2, respectively. Thus, these curves 
lack the £> = plateau. 

In addition, we note that a 8e 2 //i-step occurs at v = t 
for B = 10.1 T and B = 15.2 T. For B = 12 T, all the 
LLs are nondegenerate and hence all the steps are of the 
height of Ae 2 /h. However, if the (TV, — , +) and (o, +, -) 
LLs are degenerate at E = t, a 8e 2 //i-step appears at 
p = t. This level degeneracy happens as the magnetic 
field B = AT:t 2 /Nehv 2 F , where N = 1,2, 3, B = 10.1 T 
and B = 15.2 T fit the condition with N = 3 and N = 2, 
respectively, and thus a 8e 2 //i-step appears at p = t. We 
also find that the disappearance of a zero-Hall conductiv- 
ity plateau is always accompanied by the occurence of a 
8e 2 //i-step at p = t, because if B = irt 2 /Nehv 2 ?, B would 
satisfy B = Airt 2 / N' ehv 2 F with N' = AN. Interestingly, 
here we find that the structure of the Hall plateaus of the 
AA-stacked BLG would be significantly affected by the 
applied magnetic field, which is quite different from the 
MLG and the AB-stacked BLG. 



magnetic field is distinctly different from both the MLG 
and AB-stacked BLG as well as the conventional quan- 
tun Hall materials. In particular, the AA-stacked bilayer 
graphene could possess the unique v = plateau, in con- 
trast to other graphene materials such as monolayer and 
AB-stacked bilayer graphene. The shift of level anoma- 
lies due to interlayer hopping energy is attributed to be 
the origin of the 9 = plateau. Nonetheless, the v = 
plateau across p = would disappear if magnetic field 
B = irt 2 /Nehvp. In addition, we find that the disappear- 
ance of a zero-Hall conductivity plateau is always accom- 
panied by the occurence of a 8e 2 //i-step at p = t. Fur- 
thermore, when \p\ = [(y / rii + ^/n^) / {\Jn2 — y/rii)]t, the 
AA-stacked BLG lacks the v = ±4(m + 712)71 plateaus, 
which exist in the AB-stacked BLG. We also find that 
when r — > and H 4 U! 4 ^> (p + vt) A , a xx — >■ if the Fermi 
level is not a few Ts above and below a Dirac point. This 
implies that at the high magnetic field, the external elec- 
tric field cannot drive any current for p < t while the cur- 
rent would be perpendicular to the external electric field 
for p > t. We hope that our predicted interesting charac- 
teristics of quantum Hall effect in the AA-stacked BLG, 
which are not seen in both the MLG and BLG as well 
as the conventional quantum Hall materials, would stim- 
ulate experimental effort on fabricating the AA-stacked 
BLG and also on measurement of its transport property 
in the near future. 



IV. SUMMARY 

In conclusion, we have calculated both the quantized 
Hall conductivity and longitudinal conductivity of the 
AA-stacked bilayer graphene within linear response the- 
ory by using Kubo formula. We find that the depen- 
dence of the Hall plateau of the AA-stacked BLG on the 
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